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The Maxwell-Bloch system describes a quantum two-level medium interacting with a classical 
electromagnetic field by mediation of the the population density. This population density variation 
is a purely quantum efi'ect which is actually at the very origin of nonlinearity. The resulting nonlinear 
coupling possesses particularly interesting consequences at the resonance (when the frequency of the 
QQ ' excitation is close to the transition frequency of the two-level medium) as e.g. slow-light gap solitons 

, that result from the nonlinear instability of the evanescent wave at the boundary. As nonlinearity 

. couples the different polarizations of the electromagnetic field, the slow-light gap soliton is shown to 

^N) ' experience effective scattering whith charges in the medium, allowing it for instance to be trapped 

^ or refiected. This scattering process is understood qualitatively as being governed by a nonlinear 

Q . Schrodinger model in an external potential related to the charges (the electrostatic permanent 

' background component of the field). 

o 

C\| ■ I. INTRODUCTION 



I In the field or interaction of radiation with matter, the resonant scattering of an electromagnetic pulse with a 
p I . structure imprinted in a medium (as a periodic grating) has become widely studied for its rich and novel nonlinear 
mechanisms [IJ. Such studies have been developped in nonlinear optics where the medium is a periodic arrangement 
of dielectrics (Bragg medium) and where the nonlinearity results from Kerr effect [2], or in Bose- Einstein condensates 
where the periodic structure is produced by external applied field (thus called an optical lattice)^, or else in 
photorefractive media where an externally applied electric potential modifies locally the index value 

In a Bragg medium, the underlying model is the Maxwell equation for the electromagnetic field where the optical 
index varies along the propagation direction and where the constitutive equations (relating field and polarization) 
take into account the third order suceptibility. A simplified model is obtained then in the slowly varying envelope 
jCi^ ' approximation: the nonlinear Schrodinger equation embedded in a potential which represents the structure of the 
, medium. In a Bose-Einstein condensate, the model is the Gross-Pitaevskii equation which is again a nonlinear 
Schrodinger equation in an external potential representing the optically generated structure (the optical grating). 
' , Last, in a photorefractive medium, the variations of the applied electrostatic field induce a local variation of the 
' optical index, and possibly also of the group velocity dispersion, which eventually results again in a NLS like model 
00 I with variable coefficients. 

' It is remarkable that these situations share the same model and that it is precisely the simplest nonlinear version 
of the paradigm model for wave scattering, the Schrodinger equation of quantum mechanics, namely 



> 



+ g + HV-l^^ = (I.l) 

where the wave function ^ = il)[z^t) is submitted both to self- interaction (nonlinearity) and to scattering with the 
external potential V (depending on space z and possibly also on time t). 

In fact the common original mechanism is the coupling of an electromagnetic radiation of a given frequency cj, 
to a quantum two- level system whose transition frequency O is close to to. The interaction is then mediated by the 
variation of the population densities of the levels which makes it a nonlinear_process. The semi-classical model of 
such a process is the Maxwell-Bloch (MB) system which can be written @, 0, 



NE, (1.2) 



dN . dP 



in the case when the dipolar momentum jli2 is assumed to be parallel to the applied electric field E (at least in 
average). The constant A = (ry^ -)-2)/3 is the Lorentz local field correction, rj is the optical index of the medium and 
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the transition frequency. The dynamical variables are the real valued 3-vectors E (electric field) and P (polarization 
source), and the inversion of population density N. The above first two equations are a convenient rewriting of the 
Bloch equation of quantum mechanics for the density matrix p where one defines 

P = iVo(p2lAil2 + P12M21), ^ = ^o(P22 - Pll), (1-5) 

with A^o being the density of state of active elements (note that hcrmiticity guaranties that N and P are real valued). 
To be complete, the vector jli2 is defined from the two eigenstates \ui) of the unperturbed Hamiltonian Hq by 

Aii2 = {ui\er\u2) , H^Ui = EiUi, £2- £1 = hfl, (1.6) 

when the electric dipolar Hamiltonian is Hq ~ er ■ E. 

In the above MB system, the variations of the dynamical variable N (inversion of population density) is on the 
one side a purely quantum effect and on the other side the very mechanism of nonlinearity. Indeed, the linear limit 
(small field values) comes from assuming TV constant (in fact N ~ —Nq if the medium is in the ground state), and 
MB reduces to the Lorentz classical linear model (equation (jll.9[) below). Moreover, the structure of equation (jl.3[) 
implies effective coupling of all 3 components of the field, a property which will allow us to couple a polarized tranverse 
electromagnetic component to the longitudinal electrostatic one. 

Inclusion of a particular physical situation consists for the MB system in a convenient choice of the set of initial 
and boundary values for a particular structure of the fields (as e.g. a propagation in one given direction). The 
freedom in such choices makes the MB system a very rich model that has not finished to offer interesting results. 
Many different limits of MB have been studied, for instance the case of a weak coupling (under-dense media) allows 
for self-induced transparency (SIT) of a light pulse whose peak frequency is tuned to the resonant value fl, when a 
linear theory would predict total absorption l3| HI. The limit model in the slowly varying envelope approximation 



results to be integrable jl2l | and to have the mathematical property of transparency: any fired pulse having an area 
above a threshold evolves to a soliton, plus an asymptotically vanishing background [Ts] . Description of SIT process 
within the inverse spectral transform is fruitful here because the incident pulse, which is physically represented by a 
boundary datum, maps to an initial value problem on the infinite time line. 

Away from the resonance and for dense media, a widely used approach models the dynamical nonlinear properties 
of pulse propagation by selecting propagation in one direction. The resulting reduced Maxwell- Bloch system has again 
the nice property of being integrable [3, [l^ (also when detuning and permanent dipole are included [l^). As a 
consequence the properties of the gap soliton such as pulse reshaping, pulse slowing, pulse-pulse interactions, are 
fairly well understood, more especially as the reduced MB system possesses explicit N-soliton solutions Others 
interesting features include pulse velocity selection [11]. However, unlike the SIT model, the reduced MI3 system 
happens to be integrable for an initial pulse profile which makes it almost useless to study the scattering of an 
incident light pulse. 

To describe explicit scattering process, the reduced MB equation was replaced by the coupled-mode Maxwell- Bloch 
system where the electric field envelope contains both right-going and left-going slowly varying components (19| . 
Although the adequation of the model to the physical situation is a very difficult question often left apart, the 
approach allowed to prove of existence of gap 27r-pulses in the presence of inhomogeneous broadening, and to discover 
"optical zoomerons". Moreover numerical simulations have demonstrated the possibility of "storage of ultrashort 
optical pulses" [20]. This storage can moreover be externally managed to release the stored pulse and thus create a 
"gap soliton memory" with a two- level medium [2l| . 

The coupled-mode a ppr oach has been also ap plie d to understand the properties of resonantly absorbing Bragg 
reflectors introduced in [2^ and further studied in [23l . [23 ] , which consist in periodic arrays of dielectric films separated 
by layers of a two-level medium. Very recently, the coupled-mode MB system has been used to model "plasmonic 
Bragg gratings" (25| in nanocomposite materials where a dielectric is imbedded in a periodic structure of thin films 
made of metallic nanoparticules [265. 

Beyond such a set of approximations, the Maxwell-Bloch system offers a natural intrinsic nonlinear coupling of the 
various field components. It has been demonstrated for instance that, for a uni-directional propagation, the coupling 
of two transverse electromagnetic components can be described in the slowly varying envelope approximation by a 
(non-integrable) set of two coupled nonlinear Schrodinger equations [27i] . Thus the question of the analytic expression 
of the fundamental soliton solution is still open. However such a coupling process allows to conceive a method to 
manipulate light pulses with light by propagating a pulse on a background made of a stationary wave, when wave 
and pulse are orthogonally polarized [28,] . 

The problem of the generation of a pulse living in the stop gap of the two-level medium has been solved in [2^ by 
exciting the medium at one end with a cavity standing wave. As the frequency is in the stop gap, the result is an 
evanescent wave, at least at the linear level. It happens that the evanescent wave is nonlinearly unstable, which results 
in the nonlinear supratransmission process (soj that generates gap solitons propagating in the two-level medium at a 
fraction of light velocity. 
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We demonstrate here that the effective couphng of a tranverse electromagnetic (polarized) field to a longitudinal 
electrostatic component having a permanent background (representing local charges) allows to scatter a slow light gap 
soliton (SLGS) and to trap it in the medium or to make it move backward. This result is obtained by first studying 
the MB system in a particular situation of a propagation in a given direction (say z) of a field which is linearly 
polarized in the tranverse direction (say x) and which interacts, through the variations of the population density, with 
a longitudinal component. Then by appropriate boundary values for the tranverse component and initial data for the 
longitudinal one, we perform numerical simulations that show the propagation and scattering of the SLGS. 

We thus demonstrate that the presence of charges in a definite region of space (obtained for example by an applied 
electrostatic potential or with a doped semiconductor, or else with inclusion of metallic nanoparticles) produces a 
dynamical interaction with an electromagnetic radiation by means of a nonlinear coupling through the plasma wave 
field component (spontaneously generated out of the permanent electrostatic background). It is worth mentionning 
that this work makes use of the results of [29l w here a SLGS is generated in homogeneous two-level media and that 
the perturbative asymptotic analysis follows p7| and will not be detailed (see the appendix). Last, another way of 
manipulating light pulses have been proposed in (28l | by creating in the two-level medium a standing electromanetic 
field background (inside the passing band) which is a completely different physical process. 

The scattering is finally given a simple meaning by writing a nonlinear Schrodinger model for the envelope of the 
electromagnetic field in an external potential resulting from the electrostatic permanent background. Interestingly 
enough, the usual asymptotic perturbative analysis of the system leads to a deformed nonlinear Schrodinger system 
which allows to calculate the correct expression of the initial vacuum. However the model docs not furnish an accurate 
description of the dynamical properties of the SLGS. This is the result of the intrinsic nature of the method which 
separates the orders of perturbation when the original model naturally couples them. 



II. THE MODEL 



A. Dimensionless form. 



The system p.2IIL4[) can be written under a dimensionless form by defining first the new space-time variables 

r' = -flr, t'^nt (II. 1) 

c 

and the new field variables 



E' = ,/^E, P' = aW— P, N' = —, /ii2' = ^ (IL2) 

where by definition of the optical index e = eo'7^i a-nd where 

Wo = Nohn/2, ^2 ^ . (Ij 3) 

Note that Wq is a reference energy density. Then the MB system eventually reads (forgetting the primes) 



dt^ 



P = -aiVE, 



dN ^ dP 

_ _ ^ d^E d'^P 
Vx VxE ■ 



dt^ a<2 ■ 

We check that the unique remaining coupling constant a, defined by 

is indeed dimensionless. The length ^ of the dipolar moment is in units of C • m, the energy Ml in J while the 
permittivity e can be expressed in ■ ■ . Using finally that Nq is a density in effectively leads to a 
dimensionless a. The actual dynamical variables are the dimensionless quantities E(r, i), P(r, t) and N{f, i), namely 
7 scalar real variables, fully determined by the MB system pi.4|) with convenient initial boundary value data. Note 
that the scaled inversion of population density N now varies from iV = — 1 (fundamental) to TV = 1 (excited). 
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B. Polarized waves. 



We restrict our study to unidirectional propagation for a field which possesses a linearly polarized electromagnetic 
component E{z,t) and a longitudinal component F{z,t). Namely we assume the following structure 





E= , (II.6) 



for which the system (|II.4p reads 

Pu + P = a(l-n)E, Qtt + Q ^ a{l - n)F, 

En - E,, -Ptu Fu = -Qtu (HJ) 

nt = EPt + FQt. 

where we have defined the scaled density of excited states 

n{z,t) = 1 + N{z,t), (11.8) 

The compatibility of the chosen structure with the MB model can be demonstrated at the linear level obtained for 
iV = -1 (orn = 0) 

„ „ _ „ d'^E d^P 



Ot^ ■P = -E, VxVxE+^.-^. (II.9) 

which is nothing but the Lorentz model. The general fundamental solution possessing the structure (jll.6[) can be 
written 

^r{^t~-kz) + I I e''^"* +1 I , (11.10) 

^0 

1 c^i'^*-''^-) -I- 





l-tj2 
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The above fields are constituted of 3 basic parts: the first one is the tranverse electromagnetic component where k is 
given by the dispersion relation (plotted on fig[l|) 

c.2 = l + a, (11.11) 

cj^ — 1 

the second one is the plasma wave of frequency luq and arbitrary ^-dependent amplitude f{z), the third one is 
the electrostatic field Fo{z) that represents the permanent charges. Indeed the Gaus theorem relates Fq to the 
(dimensionless) charge density p by V • (E + P) = p = (1 + a)dzFo{z). 



C. Initial data. 



The above structure of the linear solution shows that it is possible to eliminate the variable Q{z,t) by integrating 
one equation in the system pi.7p . Actually the equation Fu — ~Qtt results from the original equation VxH = dtT) 
(no current of charges). As the propagation is along z, the third component of the curl vanishes, therefore {dtD)z = 
= dt{F + Q), which allows to express Q{z,t) as Q{z,t) = —F{z,t) + C{z), where the integration constant C{z) 
remains to be computed. It is naturally determined from the initial conditions. We consider here a situation where 
the presence of stationary charges produces a permanent electrostatic component Fq{z), namely when 

F{z,t) = (j>{z,t)+Fo{z), <j>{z,0)=0. (11.12) 

Then the structure (jll.lOp of the linear solution then implies the following expression for the longitudinal component 
of the polarization 

Q{z,t)^aFo{z)-c^{z,t), (11.13) 

where, in the absence of initial electromagnetic irradiation, we have the initial conditions {n(z,0) = 0, 0(z,O) = 0}. 
We now demonstrate that the initial vacuum {E = P = 0, n = 0, = 0} and {F — Fo(z), Q — aFo(z)} is nonlinearly 
stable as soon as Fq is below some threshold. 



k 



FIG. 1: Plot of the linear dispersion relation of Maxwell-Bloch for a — 1. The two curves are the two solutions tiJi(fc) and 
uJ2{k) of expression (lILllf) . We have also plotted the plasma frequency ujo and the free photon line lo = k. Solitons will be 
generated in the forbidden band gap ]1, ljo]. 

D. Nonlinear stability. 

The stability of the initial vacuum is proved for the dynamical variable (j){z, t) without applied external electromag- 
netic field E{z,t) (the x-component). Replacing expressions ()IL12|) and (jII.13|) in the system (jll.7|) . the equation for 
the electrostatic wave (j){z, t) becomes 

(t)tt+i^l<P^an{(t) + Fo). (11.14) 

Without applied electromagnetic wave {E = P ^ 0), the equation for n{z,t) (the density of population of the excited 
state) can be integrated taking into account n{z, 0) = and (f){z, 0) = to give 

n = -^^^-^Fo, (11.15) 

which, in the evolution for cf), provides the closed equation 

^tt + (c^o + a^o + f '^'('Z' + SFo) = 0. (11.16) 

This is an equation for an anharmonic oscillator whose frequency depends on the parameter Fo(z). The solution 
corresponding to (j)(z,0) — and (/)t(z,0) = is obviously (f>{z,t) = 0, which shows that the initial conditions are 
consistent. To study the stability, we rewrite the above oscillator as 

The condition of nonlinear stability is reached when the function £p{4>) has one single minimum in cj) = (initial 
condition). The derivative 

cj,^ + 3Fo(b + 2i^+F^) . 
a 

vanishes always at ip — 0, and there will be no other solutions as soon as maxzjfQ } < 8(1 + a) /a (in our case a = 1, 
i.e. 8(1 + a)/a — 16). With a single root at ip = 0, the potential £p{(j)) has a single minimum in (/> = 0. However, 
when max^ll^ol} exceeds the threshold value (namely the value 4 with a = 1), the potential £p{(j)) becomes a double 
well with a local maximum at point = which is thus unstable (and might be of interest to study). Thus here we 
shall consider values of Fq such that max^HFol} < 4 such that nonlinear stability condition be always satisfied. 

Thus we have obtained the initial state reached by the medium submitted to electrical potential (charges), namely 
to a permanent longitudinal component -Fo(^) which has been shown to generate a permanent polarization aFo{z) 
accordingly with (|II.13p . These data constitute then the initial condition (initial equilibrium ground state) that is 
now inserted in the model itself. 



di 
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E. Final model. 



The model studied from now on eventually reads 

Ptt + P = a{l - n)E, 
Ett — Ezz = —Ptt, 



(ptt + t^lcj) = an{cf) + Fo), 
nt = EPt-{(t) + Fo)^t, 



and it is completed by the following initial values (ground state at rest) 

F(z,0) = 0, St(z,0) = 0, P(2,0) = 0, Pt(2,0) = 0, 
(/.(^,0) = 0, 0t(z,O) = O, n(z,0) = 0. 



(11.17) 



(11.18) 



Two boundary values for the electromagnetic components are now needed, chosen following [29j as representing an 
electromagnetic cavity coupled to the two-level medium by using its stop gap as one side mirror (the cavity works at 
frequency u in the gap). We then assume an open end z — L (vanishing of the magnetic component), therefore we set 



£;(0,i) = a sin(tjt), EAL,t) = 0, 



(11.19) 



with Lu e [1, ujq]. Last, Fo{z) is a given function that results from the presence of charges, (created e.g. with a 
doped smiconductor, an applied static electric potential, inclusion of metallic nanoparticles,...), and that needs to be 
evaluated in a given physical context. We present in figOthe scheme of principle of the device that corresponds to 
our problem. 




cavity 



charges 
two-level medium 



FIG. 2: Principle of a situation where the tranverse generated electromagnetic gap soliton would scatter with a permanent 
longitudinal background Fo{z) created by a region with charges (doped smiconductor, applied static electric potential, inclusion 
of metallic nanoparticles...) where a longitudinal component F{z,t) is generated by nonlinear coupling. 



III. NUMERICAL SIMULATIONS 



Our purpose is simply here to show that a given function Fo{z) acts as a scattering potential on the dynamics of 
the slow light gap soliton. To that end one needs first to generate the SLGS. It is done by applying the principle of 
nonlinear supratransmission [sO] to the MB system as described in [2^. There it is proved that the boundary value 
pi.igp at a frequency in the stop gap, which would linearly produce an evanescent wave in the medium, is actually 
unstable above a threshold given by 



V + 3 

where lu is the driving frequency in the stop gap [1, loq]. 

Let us simply mention that the threshold of nonlinear supratransmission is obtained from the nonlinear Schrodinger 
limit of MB by assuming a driving frequency close to the gap edge uiq- For instance it is shown in (29j this threshold 
is correct to a high precision close to the edge (typically for w = 1.4 with the edge ujq = V^), precision that decreases 
to 15% at the frequency 1.2. But the very question is not the value of the threshold but its existence, source of gap 
soliton generation. Note that one needs to ensure compatibility of boundary values pi.l9|l with initial data (|II.18p 
and thus we always start the boundary driving smoothly at t > such that £^(0, 0) = 0. Note also the the stability of 
the static solution is used to determine a stable vacuum solution inside the region with charges, which is disconnected 
from the instability of the evanescent wave living in the region without charges. 
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FIG. 3: Scattering of the SLGS generated by the boundary (|II.19[) with a = 0.4 and lo — 1.4 onto the step function Fo{z) of 
(|III.2|) with zo = 50 and values of step height B as indicated. These are intensity plots of the energy density flux —EP/2 where 
the brightest regions have value 0.3. 



We display in figl3] the result of the scattering of a SLGS onto a step potential 

Fo{z) = B0{z-zo), (III.2) 

where 9 is the Heaviside function. We have used as boundary pi.lQp the following driving (at frequency lu — 1.4) 

E{0, t) 0.2[tanh(0.2(< - 20)) - tanh(0.5(i - 150))] cos{ujt), (III.3) 

The scattering mechanism will be analyzed in the next section, let us simply mention here that the SLGS is very 
robust and has been checked to (almost) maintain its amplitude and frequency across the step, the only varying 
parameter being its velocity. This property will allow us to understand the existence of a threshold step height above 
which the SLGS is reflected. Below this threshold the SLGS simply slows down to ajust to a different medium. 



IV. A PHENOMENOLOGICAL MODEL 



A. The NLS model 



To describe the gap soliton scattering in a simple and efficient way, we propose here the following model 

— 1 1 

- 2zc^o^t + V,, + -(^2 _^ 3)|^|2^ = v{z) ^. (IV.l) 

similar to p.l[) and obtained from (|A.19p by replacing the action of ip by that of an external potential V{z). The 
related initial-boundary value problem is naturally deduced from relation (IA.12P and boundary value (|II.19p as 

^{z, 0) = 0, V(0, t) = -i^e*^'^-"")*. (IV.2) 

This constitutes a phenomenological model where value of the applied potential V{z) is obtained by comparing 
numerical simulations of the above equation with those of the Maxwell-Bloch system (IIL17I) . We have obtained that 

V{z) = i H^^^o(.)^ (IV.3) 
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provides a quite accurate description of the scattering process. The above particular form has been inspired by the 
nature of the asymptotic expansion in eq. (|A.15p with the structure (|A.17p . especiaUy for the dependence in Fg , but 
the precise expression has been worked out to generate scattering processes similar to those obtained from Maxwell- 
Bloch. This is illustrated by figH] obtained by solving (jIV.ip with the above potential and boundary value. Then to 
compare to the plots of figlH we contruct the electromagnetic energy density flux —EP/2 which from (|A.12[) is given 
by 

--EP~-{il;e"^''^+^e''"^°^f. (IV.4) 
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100 



FIG. 4: Test of the soliton scattering by use of the NLS model UV. ip with boundary driving UV.S^ where a — 0.4 and u = 1.4, 
and with V{z) given by iTTOI) where Fo{z) = Be{z-AO). We have plotted the energy flux density —EP/2 given here by J/F.^l ), 
which concur with simulations of Maxwell-Bloch of fig[^ 



B. Scattering simple rules 

The interest of having such a simplified model is to derive explicit formulae for the scattered gap soliton. The 
NLS equation in a potential has been widely studied, e.g. by use of the inverse spectral theory [33], and especially 
when the potential is a Dirac delta function representing a local inhomogeneity '33] [34*1 'sB'l , which has also been done 



for discrete systems [36[. An interesting application is the study, with a stochastic scries of delta functions, of the 
competition between disorder and nonlinearity (37j [ssj . In our case the step (or barrier) height is not small compared 
to soliton amplitude (Vq is of the order of maxlipl"^) and one cannot call to perturbative approach to get analytical 
description of the scattering process. However, we shall obtain by quite elementary arguments, an expression of the 
scattered soliton velocity in terms of the incident velocity and the parameters of the medium. 

We are interested in potentials V{z) piecewise constant (in regions sufficiently larger than the soliton extension) 
for which the NLS equation pV.ip possesses approximate soliton solution in each region, far enough from the discon- 
tinuities. The soliton solution of (jIV.ip for any V{z) = Vq constant is given by 

Api{f:Sz-p,t) 

V.(z,t) = — ^- -r (IV.5) 

cosh[7(z — vt)\ 



with the following 3 relations linking the 5 coefhcients A (amplitude), /? (wave number), /i (frequency), 7 (stiffness) 
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and V (velocity). These are 



13 = ^ V, (IV.6) 



A 



7^ = ^K^ + 3)-P^, (IV.7) 
yo = (7'-/3')^^^^-2Mc^o. (IV.8) 

^0 

One important consequence of the above relations is obtained by elimination of /3 and 7 and writes 

"'^^ = 4^(^0+3) (IV.9) 

which furnishes the soliton velocity from its amplitude A and carrier frequency fj,. Note that the solution is a gap 
soliton which means /i > 0. 

Now we assume, accordingly with observations of a number of numerical simulations, that the amplitude A and 
the frequency /i do not change across a step from a region 1 where = to a region 2 where V ^ Vq. The soliton 
is produced in region 1 by the boundary driving and are measured for each numerical experiment as in p9l | (it is 
actually easier to measure v and /i and to deduce A). Then the scattering process is understood in a quite simple 
manner. First formula (|IV.9p immediately implies the existence of a particular threshold value of the constant Vq for 
which the velocity in region 2 would vanish and above which this soliton is not a solution. This threshold, called Vth, 
is thus obtained by setting v = in (jIV.9p and reads 

Vth = ^{u^l + 3) - 2cjofi. (IV. 10) 

To be clear, this threshold indicates that the NLS equation pV.l[) with Vq > Vth does not support a soliton solution 
with the given parameters A and fi. So we expect a reflection on the step of height Vq as soon as Vq > Vth- 
Below the threshold, the velocity v' in region 2 is deduced from pV.9|) and PV.10|) as 

Vo<Vth : v'^^v'-Vo"^^, Vth = v'^^. (IV.ll) 

Here the threshold value Vth given in pV.10[) is written in terms of the emitted SLGS velocity v only by computing 
hereabove the value of Vq for which v' = 0. 



C. Step potential 

In the case of the one-step potential pil.2|) we have 



and the velocity v' reads from pV.lip 

2^4 t^§ + 4 ■ 



(IV. 13) 



As a consequence the threshold amplitude Bth of the electrostatic permanent background Fq{z) — B9{z — Zq) can be 
written 

5*. = --„^V2g|±ip. (IV.14) 

We have tested these formulae against numerical simulations of the Maxwell-Bloch system (not simply the NLS 
model) , which is summarized in figure [51 Clearly the small discrepancy between soliton scattering and elementary 
formula pV.13p comes on the one side from the assumptions of conservation of shape across the step, on the other 
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FIG. 5: Soliton velocity v' beyond the step location zq in terms of the step height B according to formula f/T/. jg)) plotted as the 
continuous curve, and measured from simulations of Maxwell-Bloch (points). 

side from losses by phonon emission. Note that formula (IIV.13P relates the kinetic energy of the soliton in the two 
media through the potential energy Vq of the obstacle. Another interesting aspect of this formula is the fact that 
the scattered velocity v' depends only on the incident velocity w, once the step height Vq and the coupling constant 
a = ujq — 1 have been fixed. 

Such an expression can be confronted to numerical simulations of Maxwell-Bloch for different values of the coupling 
parameter a and different incident SLGS velocities v (obtained by varying the driving amplitude). 
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0.70 
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0.18 


0.34 


0.38 



(IV. 15) 



Some results are presented in table [W. 1 51 where a is the coupling parameter, a and to are the chosen driving amplitude 
and frequency for Maxwell-Bloch as defined in pi.l9p . Then v is the measured soliton velocity and Bth the theoretical 
threshold for soliton reflection as given by pV.14[) . The effective threshold B is then obtained by multiple simulations 
of Maxwell-Bloch where Fo{z) = B9{z — zq) varying the value of B. Such a result shows that the empirical formula 
pV.3p for the potential V{z) entering the nonlinear Schrodinger model correctly describes the scattering process 
within the chosen parameter range. 

Conclusion and comments 

The main result we wish to emphasize is the demonstration that a slow light gap soliton in a two-level medium 
can be manipulated by use of a permanent electrostatic background field. Such a property may have interesting 
applications as it is quite generic and fundamental: a two-level medium constitutes the basic model for any medium 
submitted to monochromtic radiation close to one transition frequency. 

One interesting open problem is the description of the interaction process through the method of perturbative 
asymptotic expansion. We have indeed discovered that it is the interaction of processes of different orders which is 
fundamental in the coupling of the electromagnetic field with the electrostatic component. We have not been able 
to include such interaction in the perturbative asymptotic method as by nature it decouples the different orders. 
Although the simple nonlinear Schrodinger model (jIV.l[) describes quite well the SLGS scattering of Maxwell-Bloch, 
one would appreciate to obtain it by some limit procedure. Actually the NLS equation is a natural limit, it is 
the expression of the external potential (|IV.3p representing the electrostatic permanent background that should be 
rigorously derived. 
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One may of course play with different types of potentials, as for instance a barrier of height B and width D. In 
that case at a given height, and a given incident soliton velocity, we observe tunnelling when the width value is below 
a threshold (e.g. we have obtained tunnelling with parameters of figl3]at height B = 0.8 for a width D < 4.5). This is 
an interesting problem, obviously related with the delta-function case [33|, that will be considered in future studies. 

Note finally that the soliton is generated from the boundary driving (jII.19[) by the nonlinear supratransmission 
mechanism which is an instability of the evanescent wave [soj . As a consequence, the soliton characteristics are quite 
sensitive to the driving parameters: frequency, amplitude, and shape. The figl3]has been obtained for the particular 
boundary value piI.Sp and the soliton is generated at a velocity v = 0.88. A different shape of that driving, keeping 
same amplitude and frequency, may generate a SLGS with a slightly different velocity. Note that the fundamental 
question of the relation between soliton characteristics and driving parameters is also an open problem. 



APPENDIX A: PERTURBATIVE ASYMPTOTIC ANALYSIS 




The purpose of this appendix is to provide useful! expressions for the asymptotic model derived from Maxwell-Bloch 
as in [27| but here with non-vanishing permanent fields. Although computations are lengthly, they are standard and 
we list here only the relevant results. 

a. General expressions It is actually convenient to consider the full 3- vector case and to reduce it to our particular 
interest at the very end. The dimensionless Maxwell-Bloch system (|II.4p reduced to a propagation in the z-direction 
becomes 

d^P + P - aE = -miE, 5*71 = E • dtP, 

{d!-<jd'^)E^~d^P, (A.l) 

where the singular matrix a is defined by 

f 1 o\ 

(A.2) 

The reductive perturbative expansion is a method to extract information about the variations at first order of a formal 
expansion of the fields in terms of new slow variables, see e.g. the tutorial paper [sij and the refrences therein. We 
seek slowly varying envelope approximation solutions (SVEA) for a carrier wave at frequency lo very close to loq (and 
in the gap: w < wq), thus we assume the expansion 

oo +j 

oo +j 

j = l m=-j 
oo +j 

" = E^' E N^i^,r)e^"'^"\ 

j=l m=-j 

in which the reality condition implies 

ET" = E™, p-™ = Pf, N-"'^N^. (A.4) 

An important issue is to include in the expansion the description of the permanent background Fq{z). Then we must 
assume that E and P do possess zero-frequency modes, namely that E° 7^ and P° ^ 0. We shall find that these 
quantities must obey some general constraint but are not governed by an equation (this is the reason why in general 
they are assumed to vanish) which will allow us to fix them as external data (scattering potential). The slow variables 
^ and T are here 

^ = ez, T = eh. (A.5) 



Indeed, the natural space variable is usually ^ = e{z — vt) where v is the group velocity at the frequency of the carrier. 
Here this frequency is ujq for which the group velocity vanishes, as implied by the dispersion relation pi. lip . Note 
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that this property ensures that the boundary- value problem in the physical space maps effectively to a boundary value 
problem in the slow space. 

Inserting the infinite series expansion (jA.3|l into the system (jA.ip we obtain, after careful (tedious) algebraic 
computations, the following closed form system at order (remember Ej"^ = E}) 



2i'^drEl + ^a|E} = -iV°E} + (E} • E;)E; + ^(E} • E})E-i, 



a 



a 

UJ, 



1 



a.iv°-a.(fE;.E? + i%±|E}.Er^ 



(A.6) 



together with the following constraints on the backgroung field E^ 

/ a| 

9| I E'l^ = 0. 

V 92 



(A.7) 



The principle of the method is then to express all relevant quantities in terms of the solution of the nonlinear evolution 
(IA.6p . These are given by 



E =e (E° + Eje*"°* + E^^e*'^''*) + ©(e^) 
P =e (aE? - E}e"^«* - Ej^^e"*""*) + 0{e^) 

1. 
2" 



n ^e^{N^ - E} • E°e*"°* - -E} • Eje^*""* - E^^ • Eje" 



■iiijQt 



(A.8) 
(A.9) 



(A.IO) 



where E} and iV^ are given in equations (|A.6|) . The structure of the external datum Ej resulting from the constraint 
(|A.7P reads 



E? = 



(eUt) 



(A.ll) 



where the dependences can be arbitrarily fixed. 

b. Electrostatic potential By inverse scaling we may come back to the original physical variables z and t, and 
we restrict our study to the case (jll.6l) where the longitudinal component F(z,t) = 4>{z,t) + Fq{z). Rewriting then 
expressions (jA.8|) - (IA.10p in physical variables, we get at first order (namely e for E and P, for n) 



E 












+ 







e'"^"* + c.c. 














( H^^,t)\ 











e'"^"* + c.c. 






\^{z,t)} 





Fo(pe*""'* + -(V'^ + Lp^) e^*"'"* + c.c. 
where we rename {E\)x = ip, (El)z = (f and TVj = The corresponding evolutions (|A.6|) eventually read 



(A.12) 

(A.13) 
(A. 14) 



'2iLUoipt H ^ — iJJzz = -avijj + — {-tj/ + (p^)ii, 



^0 + 1 



(A.15) 
(A.16) 
(A.17) 
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Thanks to the initial data (jlLlSp . namely here to -0(2, 0) = 0, 0) = and 0) = 0, we may integrate now the 
evolution of v to get 

-=4^(H^ + M^). (A.18) 
Inserted in the equations for -0 and 93 it gives the final NLS-like system 

-2ic^o^t+^^^^^ + ^(^^0 + 3)|i^| V = 

iK2-iV2^-(,.2^1)|^|V, (A.19) 

iK2_i)^2-_(^2^^)|^|2^_ (A.20) 

The equation (jA.19p is a nonlinear Schrodinger equation for i/j (envelope of the tranverse electromagnetic component) 
in some external •potential created by the envelope of the plasma wave. The applied permanent background Fq(z) 
does not act directly on ip but through the plasma wave (p according to the dynamical equation (jA.20p . 

The above asymptotic analysis has allowed to demonstrate that the elimination of the variable Q{z, t) by equations 
(III.12P and pi.l3p . namely by comes naturally as the first order solution. Then the initial condition 0) = n(z, 0) = 
was shown to be stable in the absence of applied electromagnetic transverse component. By means of (|A.12p . such 
a set of initial data maps to (/?(z, 0) = 0. The problem is that the evolution (jA.20p has the unique solution if{z, t) = 
for initial data 93(2, 0) = and ^pt{z^ 0) = (representing an initial equilibrium state). Thus wc are not able to use the 
system of NLS-like equations (|A.19[) (|A.20[) to describe the gap soliton scattering with purely vanishing initial data. 

This comes from the fact that the scattering process in Maxwell-Bloch is initiated at higher orders by nonlinear 
coupling. Indeed, considering the MB system (III.17|) with initial-boundary value problem (jll.lSP (III.19|) , it is clear that 
the interaction of E (electromagnetic component) with the permanent electrostatic background Fa{z) works with the 
birth of the plasma wave 4>. And grows on an initial vacuum by mediation of the normalized density n of population 
of the excited level, which is an effect of second order (see the time evolution of n{z, t)). But the method of multiscale 
asymptotic series separates the orders as indeed, at each order e" an equation is obtained (independent of e). We thus 
cannot expect the method to describe correctly such a situation. 

We have pursued the series to order to check that indeed no source term appear. As expected from the method, 
we got a linear system for the next order expansion with the variable coefhcients (p and 0, given from the preceding 
order. This confirms that the method does not allow feedback interaction of higher orders to the fundamental one 

Aknowledgements. This work has been done as part of the programm GDR 3073 Nonlinear Photonics and 
Micro structured Media. 
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